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Abstract 

We present resistivity measurements of Mott insulator NiS 2 single crystals after heat treatment. We find a strong increase of the 
low temperature resistivity that relaxes back towards the pristine behaviour over several days with a time constant of 45 h at room 
temperature. The low temperature resistivity has previously been shown to be dominated by surface conduction [T]. Consequently, 
the changes induced by heat treatment are attributed to changes to surface states. Our results suggest the creation of vacancies in 
the surface that re refilled from the bulk via atomic diffusion. We estimate a diffusion constant of D =» 1 x 10 '®m^ s * at room 
temperature. We identify sulphur vacancies as the most likely to form via oxidisation of sulphur forming volatile SO 2 during heat 
treatment. Our results point towards these sulphur vacancies to be the source of surface state localisation in NiS 2 . 
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1. Introduction 

NiS 2 is one of the best examples of a Mott-Hubbard metal- 
insulator system EISlIll. It can be understood analogous to 
NiO, the canonical Mott-Hubbard insulator Q: The sulphur 
atoms in NiS 2 form dimers with a total valence of two 
(cf. Fig. 0. As a consequence, Ni is divalent with a half-filled 
Cg band as it is in NiO. It is these Ni states that are split by corre¬ 
lations into the lower Hubbard band « 1 eV below the chemical 
potential n and the upper Hubbard band « 3.5 eV above fi El. 
Whilst the sulphur bonding bands pa- are well below the lower 
Hubbard band the antibonding band p* is centred at ^ 1 eV 
above p. Thus the charge gap forms between the lower Hub¬ 
bard band and the sulphur antibonding band classifying NiS 2 
as a charge-transfer insulator Q. 

Magnetic order is found at low temperatures in NiS 2 : Whilst 
NiS 2 is paramagnetic at room temperature, it orders antifer- 
romagnetically below around 38 K and eventually becomes a 
weak ferromagnet below 30K jSl. 

The bulk insulating state in NiS 2 can be tuned by either Se 
substitution or pressure. This presents NiS 2 as a prototypi¬ 
cal system to study Mott-Hubbard physics where fundamental 
questions can be addressed. In fact, recent quantum oscillation 
measurements at high pressures succeeded in demonstrating the 
emergence of a large Fermi surface encompassing the full elec¬ 
tron count in the metallic high-pressure phase a. 

The insulating state in bulk NiS 2 at ambient pressure has a 
gap of Ab ~ lOOmeV as seen in optical measurements and 
in activated behaviour in resistivity measurements at high tem¬ 
peratures Goiini. This gap is derived from the above band 
positions by taking into account the band width of the lower 
Hubbard band and the sulphur antibonding band. Surprisingly, 
a deviation from activated behaviour is seen below 200 K; the 
resistivity p{T) shows a shoulder at 150 K with a plateau at 
^ 100 K and a second strong increase of p(T) below 20 K. 



It has been suggested that the complex temperature depen¬ 
dence of the resistivity originates from a surface layer with 
metallic behaviour above 50 K and insulating behaviour below 
50 K in addition to the activated behaviour for bulk conduction 
m. Hall effect measurements show a diminishing charge car¬ 
rier concentration as the temperature is lowered from room tem¬ 
perature to lOOK in accordance with activated bulk behaviour. 
A saturation of carrier concentration below 100 K is attributed 
to surface conduction dominating at low temperatures. Sur¬ 
face conduction has been underpinned by comparison of sam¬ 
ples with different ratios of surface area to volume showing the 
expected shift in the respective contributions m. The surface 
carriers have a rather low mobility of a; 1 cm^ s“' but a 
large surface charge carrier concentration ns ~ 5 x lO''^ cm“^ 
amounting to =» 1 electron per unit cell at the surface. This in¬ 
dicates ai 1 broken bond per unit cell to generate the amount 
of surface charge carriers observed. The steep increase in re¬ 
sistivity at lowest temperatures T < 30 K has been attributed 
to Anderson localisation due to the increased effect of disorder 
m. The precise origin of the charge carriers on the surface 
has not been identified previously. 
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Here, we report resistivity measurements in combination 
with heat treatments to our sample. We observe an enhance¬ 
ment of the low temperature resistivity. This enhancement re¬ 
laxes back towards the behaviour of the pristine sample over 
several days. These results point towards a sulphur depletion 
during the heat treatment that is replenished from the bulk 
through atomic diffusion. Our results suggest that sulphur 
atoms provide charge carriers to the conducting surface state 
whilst sulphur vacancies act as trapping potentials of « 1 meV. 
Furthermore, our results suggest that it is these vacancies that 
cause the localisation of surface states at lowest temperatures. 


2. Experimental Details 


Single crystals were grown in two steps following the pro¬ 
cedure in IHJII : First, polycrystalline NiS 2 was produced via 
solid state reaction from stoichiometric ratios of high purity 
nickel and sulphur. Single crystals were then obtained from 
this starting material using the Tellurium flux technique with 
slight excess of sulphur to compensate for losses due to the high 
vapour pressure of sulphur. The resulting millimeter sized crys¬ 
tals were characterised with powder x-ray diffraction to have a 
lattice constant of 5.6893(5) A. Structure refinement of single 
crystal diffraction measurements confirmed full occupancy of 
sulphur and nickel sites in the bulk Q. 

Resistivity measurements were performed in a ^He-cryostat 
using 4-point low-frequency lock-in technique with a constant 
current of lOOnA which was tested not to cause Joule heat¬ 
ing. Contacts were made using silver loaded epoxy (Epo-Tek 
EE129-4). 

Heat treatment of the sample was carried out at 150 °C for 
2 h in a constant flow of « 0.5 ml s ' of Argon gas with a purity 
of 99.998 %. However, even this high level of purity can pro¬ 
vide enough oxygen to oxidise the sulphur of the surface of our 
crystal as estimated in Appendix A| 


3. Results 

Our resistivity measurements of the pristine sample (black 
curve in Eig.|^ show typical behaviour for NiS 2 : The resistivity 
increases towards lower temperatures with activated behaviour 
between 200 K and 300 K and with a corresponding activation 
energy of A = 104meV in good agreement with earlier reports 
of high-temperature behaviour uniiiiiiKni. Deviation from 
the high-temperature activated behaviour is attributed to surface 
conduction (cf. orange line in Eig.|^ HI: The shoulder in p{T) 
around 150 K and the plateau between 50 K to 100 K arises due 
to metallic behaviour of the surface whilst the surface conduc¬ 
tion turns insulating at low temperatures leading to the strong 
increase in resistivity below 50 K. 

The surface to volume ratio of our sample is significant for 
the rectangular cross section of 0.5 mm x 0.2 mm. This is 
in accordance with a sizeable surface contribution associated 
with the large deviation from the high-temperature activated 
behaviour below 200 K. Eor the pristine sample we observe 
metallic behaviour of the total resistivity dp/dT > 0 around 
85K(cf.Eig.[g. 



Figure 2: Temperature dependence of the resistivity of NiS 2 for our sample 
before and after heat treatment at the times given in the legend. These times 
represent time at room temperature in accordance with our model of atomic 
diffusion. Curved arrow indicates the time evolution after heat treatment. Inset 
shows same data on a logarithmic resistivity scale versus reciprocal temperature 
scale. Orange lines in inset and main panel represent the same fit to activated 
behaviour p oc exp(A/^B7') above 200 K. Straight aiTow in the main panel 
indicates the transition temperature of the ferromagnetic order. 


A peak is observed in p{T) at Tc = 30K marking the tran¬ 
sition into the ferromagnetic state E). No clear signature is 
seen in p{T) for the antiferromagnetic transition. Only in the 
derivative dp/dT a peak is visible at Tn = 38K (Eig.l^ which 
we assign to the antiferromagnetic transition ms El II. 

We now turn to our results after heat treatment in argon for 
2h at 150 °C. The resistivity before and after heat treatment 
are compared in Eig|^ Immediately after heat treatment the re¬ 
sistivity reveals a large enhancement of more than an order of 
magnitude at low temperatures whilst the high-temperature be¬ 
haviour above 200 K remains the same. This low-temperature 
enhancement is found to decrease over the course of several 
hours to days. Such a slow relaxation is a first indication of 
atomic diffusion as an underlying cause. Atomic diffusion is 
exponentially diminished at low temperatures; thus, we count 
time the sample spent at room temperature between subsequent 
resistivity measurements presented in Eig.|^ In fact, when stor¬ 
ing the sample at low temperatures {T 9 K) immediately af¬ 
ter the heat treatment the resistivity remained constant to better 
than 0.5 % over 3 h compared to a change of almost 50 % of the 
low temperature resistance after 0.5 h storage at room tempera¬ 
ture. This corresponds to a reduction of diffusion of more than 
100 fold. 

The low temperature insulating behaviour shows a qualita¬ 
tive change after heat treatment. In general it follows activated 
behaviour but displays sample dependences. Many samples can 
be fitted with an exponential form including a power-law form 
of the temperature in the exponential term 

Such a general form is predicted for variable range hopping 
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Figure 3: Derivative of the resistivity dp/dF as a function of temperature around 
the magnetic phase transitions. Data after heat treatment are offset for clarity 
with zero for each dataset indicated on the right. Arrows indicate peaks at¬ 
tributed to the onset of weak ferromagnetism at Tq and antifen'omagnetism at 

T’n. 



Figure 4: Low temperature behaviour of the resistivity. Coloured lines repre¬ 
sent fits of eq. [T]to the data between 10 K and 20 K. Data are plotted with the 
resistivity on a logarithmic scale and the temperature on a reciprocal scale. On 
this scale an exponent a: = 1 yields a straight line whilst x < \ and jc > 1 
manifest as downward and upward curvature, respectively. 


m. For our sample the pristine resistivity can be fitted with 
eq. [T]below 20 K where we find an exponent x = 0.8 and a gap 
As = 1 meV. This exponent deviates from the value of x = 1/2 
originally predicted for 2D systems M- 

Heat treatment is found to alter the low temperature be¬ 
haviour of the resistivity. Whilst p(T) can still be fitted using 
(eq. 0 the exponent becomes larger and does exceed x > 1 as 
can be seen from the upwards curvature towards lower temper¬ 
atures in Fig. compared to the slight downward curvature of 
the pristine sample. 

Strong Coulomb repulsion is expected to enhance the expo¬ 
nent in the variable range hopping model in 2D beyond x - 1/2 
but is expected to be limited to x < 1 lfT9l l20l 1211 . Thus the ex¬ 
ponent found for the pristine sample can be understood within 
the correlation corrections of the variable range hopping model 
but the behaviour after heat treatment is incompatible. This de¬ 
viation remains elusive. 

Given the complex behaviour of the resistivity with both bulk 
and surface contribution, and the surface behaviour only par¬ 
tially understood, we aim to analyse the time dependence of the 
resistivity without restricting to a specific form p{T). We take a 
Drude model as the starting point with the conductivity 

ne^ 

■ ( 2 ) 

mv 

Here, both the charge carrier concentration n and the scattering 
rate v may be temperature dependent for instance with activated 
behaviour n oc exp(A/kBr) and power law temperature depen¬ 
dence V <x. TP, respectively while e and m denote the electron 
charge and mass, respectively. 

The enhancement of the resistivity after heat treatment at low 
temperatures roughly following activated behaviour (eq. 0 sug¬ 
gests a reduction of the charge carrier concentration to be dom¬ 
inant. Two main scenarios are considered which are analysed 


with the same procedure; 

i) The reduction of surface charge carrier concentration can 
arise from trapping of available states in impurity potentials of 
magnitude Ai around vacancies created by heat treatment. This 
leads to a reduced surface charge carrier concentration 



where the vacancy concentration ny decays over time t with a 
relaxation time r. The surface charge carrier concentration is 
subject to activated behaviour with a gap As as derived above 
leading to a combined gap of the impurity and surface conduc¬ 
tion 


n — no exp 



- My exp 



Ai + As 
ksT 


(4) 


ii) Alternatively, the charge carrier concentration can be di¬ 
rectly reduced. Surface states in NiS 2 are suggested to arise 
from dangling bonds of individual atoms. Heat treatment re¬ 
moving these atoms can lead to a reduction of surface charge 
carrier concentration 



(5) 


where the surface charge carrier concentration is subject to the 
surface activated behaviour with gap As. Eq. |^can be rear¬ 
ranged to a similar form to eq. 

"="»“p “p (■;) “p (■ ^) ■ ® 


Substituting the time dependent charge carrier concentration 
of either (eq.[^ or eq.j^into the Drude conductivity we obtain 


(cr - o-q) = - 


e^ny 

mv 



(7) 
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Figure 5: Relaxation of the conductivity: The difference \cr - (TqI of the con¬ 
ductivity with respect to the initial conductivity cro is plotted against the time t 
the sample spend at room temperature after heat treatment. Solid lines are lin¬ 
ear fits with intercept a and slope b to the data in the logarithmic representation 
of the conductivity difference. 


where the gap is A = Ai + As or A = As for scenario i) or ii), 
respectively. In both cases the initial conductivity ctq is associ¬ 
ated with the initial charge carrier concentration no before heat 
treatment. This allows us to analyse the data for the two sce¬ 
narios at the same time where we will identify indications for 
scenario ii). Additional (bulk) contributions to the conductivity 
eq. [7] can be absorbed into co in a two-band model of surface 
and bulk conductivity. 

The exponential nature of the decay and gap are best anal¬ 
ysed from a logarithmic representation 


ln|cr- cTol = In 


e^riy \ 
mv I 


A t 
ksT T 


( 8 ) 


Such a logarithmic (base 10) representation is shown in Fig. 
where the initial conductivity cro and time dependent conduc¬ 
tivity cr is obtained as the inverse of the resistivity of the pristine 
sample and after time t the sample spent at room temperature 
subsequent to heat treatment, respectively. 

We find the data to be well described by this logarithmic form 
as shown with solid lines in Fig.|^in agreement with our model 
of relaxing impurity concentration. The intercept a and slope b 
in Fig. [^represent the first two and last term in eq. respec¬ 
tively: 


a - ln( ] 

intercept 

(9) 

\ mv 1 kbI 

slope 

(10) 


T 


We plot a and -1/b in Fig. as a function of temperature. Re¬ 
sults in Fig . have been converted to a natural base e in line 
with the natural logarithmic form in eqj^ Thus they represent 
a and -1/b as associated in eqs.[^and[T^ 

The temperature dependence of a and b in Fig.|^(a) and (b) 
display distinct temperature dependencies that can be linked to 
the behaviour of |(To - cr| in Fig. 0 and p(T, t) in Fig. The 
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Figure 6: Temperature dependence of the intercept a and inverse (negative) 
slope -Ijb from fits of ln|f7 - (Tq\ as shown in Fig.Solid red line in (a) 
represents a reciprocal divergence a{T) = a^- {hjk^T) between 9 K and 20 K. 
Solid line in (b) marks the constant relaxation time observed above 50 K. 
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intercept a reflecting the difference |(To - cr| immediately after 
heat treatment has large negative values at lowest temperature 
corresponding to a large difference in cr-cro and conversely the 
large enhancement of the resistivity p-po as seen in Fig. 12 The 
intercept a increases reaching a maximum at 50 K above which 
it decreases again. These two regimes are naturally associated 
with the two regimes of conduction in the surface: insulating 
below 50 K and metallic above 50 K. 

Below 50 K in the insulating regime the second term in eq. |2 
is dominant governing the trapping of charge carriers. In fact, 
we find a(T) to be compatible with a negative inverse tempera¬ 
ture dependence as derived for the second term in eq.l^ A fit of 
an inverse form yields A 1 meV as shown in Fig.l^a). This 
value agrees with the activation energy As found for the surface 
conduction in the insulating regime within our accuracy. Thus, 
this points towards scenario ii) being realised where A = As. 

The value of a reaches a maximum around 50 K and de¬ 
creases at higher temperatures. Whilst a non-monotonic be¬ 
haviour can be understood from a temperature dependent scat¬ 
tering rate v entering eq. |2 the main strong downturn above 
100 K is a signature of the diminishing surface conductivity in 
the metallic regime. In addition, the total conductivity is dom¬ 
inated by the bulk at higher temperatures resulting in the en¬ 
hanced uncertainty above 150 K. 

We now turn to the slope b representing the inverse relax¬ 
ation time as extracted from fits to the relaxation of the con¬ 
ductivity. The relaxation time is expected to be independent 
of temperature in eq. Indeed a constant relaxation time 
of T = 45 h is found above 50 K. Below 50 K we find an in¬ 
crease towards lower temperature which cannot be understood 
within our simple model. This is the regime where the sur¬ 
face contribution shows insulating behaviour. A variable range 
hopping model may be necessary to describe the behaviour in 
this range. Flere, the time dependent impurity concentration 
alters the hopping distance and possibly also the energy spec¬ 
trum of the hopping sites. Additionally, the presence of strong 
Coulomb correlations in NiS 2 may need to be taken into ac¬ 
count as demonstrated by the unconventional exponent x in the 
activated low-temperature regime. 

4. Discussion 

We model the relaxing behaviour of the resistivity with a re¬ 
duced charge carrier concentration due to vacancies created by 
the heat treatment in the surface layer. Pyrites like NiS 2 are 
well known to generate SO 2 when heated (firing pyrites is an 
established industrial procedure to produce SO 2 ). Thus, sul¬ 
phur vacancies are the most likely impurities created by heat 
treatment. 

We find the activation energy A for the missing charge carri¬ 
ers associated with the vacancy concentration ny in eq.jTjto be 
the same as the activation energy As for surface conduction at 
lowest temperatures. This points towards a low level of vacan¬ 
cies to be present in pristine samples of NiS 2 which could arise 
during the growth or by storing crystals at room temperature. 
In this case measurements of low-temperature cleaved samples 


could allow detection of surface conducting in the absence of 
vacancies on a stoichiometric sample. 

The surface conduction is dramatically reduced at lowest 
temperatures. From this we conclude that the impurity con¬ 
centration must be of similar magnitude as the surface carrier 
concentration which was estimated to be of the order of one per 
unit cell of NiS 2 on the surface. Assuming an impurity concen¬ 
tration of ny ~ 1x10^^ cm“^ in the surface layer of unit cells 
we can estimate the diffusion coefficient for the sulphur to dif¬ 
fuse from the bulk to the surface, where we assume sulphur to 
diffuse from « 10 unit cells deep. We find a diffusion coefficient 
of Z) ^ lx 10 '°m^ s' which is large but not unprecedented 
with a similar diffusion coefficient of 8 x 10" m^ s ' for Cu in 
Si at room temperature. An Arrehnius analysis of the diffusion 
coefficient as a function of temperature the sample is kept at 
can reveal the activation energy and the mobility of the sulphur 
atoms. This Arrehnius analysis may be carried out in the pro¬ 
tected atmosphere of the cryostat and below room temperature 
in order to avoid additional sulphur depletion. 

5. Conclusion 

In summary, we have discovered an enhancement in the re¬ 
sistivity of NiS 2 induced by heat treatment. This enhancement 
is found to decay with the resistivity relaxing back towards the 
pristine behaviour with a time constant of t = 45 h at room tem¬ 
perature. We identify the changes to be restricted to the surface 
conduction of NiS 2 present at low temperatures. Our analysis 
suggests sulphur vacancies forming due to oxidisation of sul¬ 
phur during heat treatment. These vacancies are refilled from 
the bulk through atomic diffusion. We estimate a diffusion co¬ 
efficient D ai 1 X lO '^m^s '. This diffusion is expected to 
follow Arrehnius behaviour which can be tested via the depen¬ 
dence on temperature at which the sample is stored between 
low-temperature measurements. 

Our results point towards sulphur atoms to provide surface 
charge carriers in NiS 2 as we find the same activation energy 
for the pristine sample at low temperatures and for the restored 
charge carriers after sulphur diffusion. This suggests that acti¬ 
vated behaviour of the surface conductivity at lowest tempera¬ 
tures is associated with sulphur vacancies present on the surface 
of pristine NiS 2 . 

The time dependent vacancy concentration on the surface 
provides a route to study variable range hopping on the surface 
of a correlated insulator where the length scale of the hopping 
can be tuned by heat treatment on a single sample. 

In addition, our results are relevant to understand and de¬ 
sign the properties of nano-crystalline NiS 2 and pyrite materials 
currently explored for applications in capacitors and batteries 

Author information and Data management Correspon¬ 
dence and requests for materials should be addressed to S.F 
(Sven.Friedemann@bristol.ac.uk). Raw data will be made 
available from the Data repository at the University of Bristol 
1221 . 


5 


Acknowledgements The authors would like to thank Chris Bell 
for fruitful discussions. We acknowledge partial support by the 
EPSRC. 


Appendix A. Estimation of oxygen present during heat 
treatment 

We performed the heat treatment of NiS 2 in high purity at¬ 
mosphere in an effort to reduce the effect of oxidisation. Yet, 
we estimate the amount of residual oxygen to be sufficient for 
full oxidisation of sulphur on the surface of NiS 2 . After ini¬ 
tial flushing the furnace we used a flow of 2 ~ 0.5niLs“^ 
This amounts to a total volume of 41 passed via the sam¬ 
ple during the heat treatment. We used 99.998% pure argon. 
Assuming the residual impurity to be oxygen we obtain an up¬ 
per limit for the amount of oxygen present to the sample of 
N 02 ~ 4 X 10 ® mol. The sample has a surface area of less then 
1 mm^. This is equivalent to ^ 1 x 10 ** mol of NiS 2 on the sur¬ 
face. This amount of sulphur is 5 orders of magnitude smaller 
than the amount of oxygen present in the used gas. Thus, we 
conclude that oxygen levels are more than sufficient to cater for 
a surface reaction in our sample. 
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